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Abstract 

The Bardeen model describes a regular space-time, i.e. a singularity-free black hole space-time. 
In this paper, by analyzing the behavior of the effective potential for the particles and photons, we 
investigate the time-like and null geodesic structures in the space-time of Bardeen model. At the 
same time, all kinds of orbits, which are allowed according to the energy level corresponding to the 
effective potentials, are numerically simulated in detail. We find many-world bound orbits, two- 
world escape orbits and escape orbits in this spacetime. We also find that bound orbits precession 
directions are opposite and their precession velocities are different, the inner bound orbits shift 
along counter-clockwise with high velocity while the exterior bound orbits shift along clockwise 
with low velocity. 
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I. INTRODUCTION 



The Bardeen model describes a regular space-time using the energy-momentum 

tensor of nonlinear electrodynamics as the source of the field equations and it is also known as 
a regular black-hole solution which obeys the weak energy condition. This global regularity of 
black hole solutions is quite important to understand the final state of gravitational collapse 
of initially regular configurations. When ratio of mass to charge is 27g 2 ^ 16m 2 , the Bardeen 
model represents a black hole and a singularity-free structure [3|]. When 27 g 2 = 16m 2 , the 
horizons shrink into a single one, which corresponds to an extreme black hole such as the 
extreme Reissner-Nordstrom solution. The physically reasonable source for regular black 
hole solution to Einstein equations has been reported around 1998 #tZ|. In the Bardeen 
model, the parameter g representing the magnetic charge of the nonlinear self-gravitating 
monopole[3|, was studied later on. 

It is well known that many effects, such as bending of light, gravitational time-delay, 
gravitational red-shift and precession of planetary orbits, were predicted by General Rela- 
tivity. Because these gravitational effects are very important for theories and observations, 
many theoretical physics and astrophysics are interested in investigating them for different 
gravitational systems. The geodesic structure with a positive cosmological constant was in- 
vestigated by Jaklitsch et al.|8j, the corresponding effective potential was analyzed in detail. 
The analysis of the effective potential for null geodesies in the Reissner-Nordstrom-de Sitter 
and Kerr-de Sitter space-time was carried out in Refs. [9]] and fll |. All possible geodesic 
motions in the extreme Schwarzschild-de Sitter space-time were investigated by Podolsky 



10j. Lake investigated light deflection in the Schwarzschild-de Sitter space-time (12J. Ex- 



act solutions in closed analytic form for the geodesic motion in the Kottler space-time 
were considered by Kraniotis et al [ijj]. Kraniotis 15] investigated the geodesic motion of 
a massive particle in the Kerr and Kerr (anti)de Sitter gravitational field by solving the 
Hamilton acobi partial differential equation. Cruz et al.(l6j studied the geodesic structure 



of the Schwarzschild anti-de Sitter black hole. Chen and Wang jl3|, Il7l-ll9| have investigated 
the orbital dynamics of a test particle in gravitational fields with an electric dipole and a 
mass quadrupole, and in the extreme Reissner-Nordstrom black hole spacetime. The mo- 
tion of test particle in Horava-Lifshitz black hole space-times was studied using numerical 
techniques [20]. 
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To find all of the possible orbits which are allowed by the energy levels for time-like and 
null geodesic in Bardeen spacetime, we analysis the effective potentials in detail. To describe 
the trajectories of massive and null particles, we have a direct visualization of the allowed 
motions. This paper is organized as follows: In Section II, we give a brief review on the 
Bardeen spacetime. In Section III, we give out the motion equations, and define the effective 
potential. In Section IV and V, we discuss the time-like and null geodesic structure of the 
Bardeen spacetime in detail. A conclusion is given in the last section. 



II. THE BARDEEN SPACETIME 



The line element representing the Bardeen spacetime is given by 
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where the parameter g represents the magnetic charge of the nonlinear self-gravitating 
monopole [3J . The corresponding lapse function is 



fir) = 1 - 



2mr 2 



(r 2 + g 2 ) 2 

The Bardeen model describes a regular space-time for the following inequality: 

2 ^ 16 2 
q ^ — m . 

y 27 



(2) 



(3) 



When g 2 < ^rn 2 , there are two horizons in Bardeen spacetime. For the equality g 2 = ^m 2 , 
the horizons shrink into a single one, which are showed in Fig.l in detail. 




FIG. 1: Horizons of the Bardeen spacetime. 
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III. GEODESICS EQUATION 



It is well known that the Euler-Lagrange equations for the variational problem associated 
to spacetime metric describes the geodesies. So we set up the corresponding Lagrangian 
according to Eq.flT]) 

C = -f(r)i 2 + firy 1 ?- 2 + r 2 (9 2 + sm 2 #0 2 ), (4) 

in which the dots denote the derivative with respect to the affine parameter r. The Hamil- 
tonian motion equations are 

n,-f = o, (5 ) 

where U q = dC/dq is the momentum to coordinate q. Since the Lagrangian is independent 
of (t,4>), the corresponding conjugate momentums are conserved, therefore 
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U t = -{!-— -r)t = -E, (6) {eq:E} 

[r z + g A ) 2 

11^ = r 2 sin 2 9<j) = L, (7) {eq:L} 

where E and L are motion constants. 
From the motion equation for 9 

n„-| = o, (8) 



we obtain 



d(r 2 6) 



r sindcosOcj) . (9) 



dr 

If we simplify the above equation by choosing the initial conditions 6 = tt/2, 6 = and 
9 = 0, the Eq.flTJ) becomes 

11^ = r 2 <j) = L, (10) {eq:Ll} 

from Eqs-OHl [7]), the Lagrangian (J3J) can be written in the following form 

E 2 f 2 L 2 
2£ = h = - ^---j. (11) 



5" ^ — T 

(r 2 +g 2 )? (r 2 +g 2 )? 



Now we solve the above equation for r 2 in order to obtain the radial equation, which allows 
us to characterize possible moments of test particles and explicit solutions of the motion 
equation of test particles in the invariant plane 

r 2 = £ 2 -(l ^^)(h+^), (12) {eq:motion} 



(r 2 + g 2 ) 1 -- 



r 



It is useful to rewrite the above motion equation as a one-dimensional problem 



• 2 



E - V. 



effi 



where V 2 ^ is defined as an effective potential 

2mr 2 



V 2 ff 

eff 



1 - 



(r 2 + g 2 ) 2 



(13) 



(14) {eq:effecti 



IV. TIME-LIKE GEODESIC STRUCTURE 



For time-like geodesic h = 1, the corresponding effective potential becomes 



V. 



eff 



2mr 2 w L 2 . 
[r z + g z )2 r 



and the orbit equation for massive particle is 





(15) 



FIG. 2: The behavior of the effective potential of non-radial particle for fixed L = 3.5m(left) and 
fixed g = 0.87m(right). 



r 2 = E 2 — (I 
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(16) {eq:time-li 



By using Eq.( flO|) and making the change of variable u 1 = r, we can obtain orbit equation 
for massive particle 

du. 2 E 2 — 1 2 2mu + 2mu 3 L 2 

-u H — _ a , (17) {eq:time-li 



d<f>' L 2 L 2 (l + ttV, 

Differentiating (flTl) . we have its second order motion equation 

2m 



3 J 
■2 



d 2 u 3mu 2 L 2 + 3m 

-T7T7 + U 
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L 2 (1 + mV)5 L 2 (1 + mV 
We solved (JT7J) and fll8p numerically to find all types of geodesies and examine how the 
parameters influence on the timelike geodesies in the space-time of Bardeen model in detail. 

From the effective potential curve (see Fig.3), we can identify 3 classes of orbits: i.e. 
planetary orbits, escape orbits and circular orbits when the energy of particle E satisfies 
two critical values E Cl and Eq 2 - 



(18) {eq:time-li 




FIG. 3: The behavior of the effective potential in Bardeen space-time with g = 0.70m, L = 3.5m, 
m = 1 and energy levels = 0.90 and Eq ± = 0.92. 

A. Time-like bound geodesies 




FIG. 4: The behavior of the effective potential of time-like bound geodesies for E 2 = 0.91,5 = 
0.7m, L = 3.5m, m = 1 

In Fig. 4 the dashed line denotes the value of the energy E 2 = 0.91, i.e. < E 2 < Eq . 
From potential curve, we can find two kinds of bound orbits for this energy level: 

I) The particle orbits on a many- world bound orbit between the range < r < Tb, 
which is near the singularity and can cross the two event horizons. The ta and tb are 
the perihelion and aphelion distance of the planetary orbits, respectively. We also can find 
the clockwise precession of planetary orbits which it is a well-known gravitational effect in 
general relativity theory. 

II) The particle orbit is on a two- world bound orbit in the range rjj < r < r^, where 
the td and te are the perihelion and aphelion distance, which are larger than the orbit of 
Case I. The two-world bound orbit is outside the event horizon. However, we also can find 
that the precession direction of the planetary orbit is counter-clockwise, and the precession 
velocity is slower than the orbit of Case I. 

These two kinds of two-world bound orbits are simulated in Fig. 6. 
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FIG. 5: Examples of the two-world bound orbit in the Bardeen space-time with E 2 = 0.91, 
g = 0.7m, L = 3.5m and m = 1. 

B. Time-like circle geodesies 

From Fig. 6 and 7, we can see that there are two different circular orbits. One is a unstable 
circular orbit, the other one is a stable circular orbit. 

I) When the energy of particle E equals the peak value Ec 2 of the effective potential, the 
particle can orbit on a unstable circular orbit at r = rc 2 . Any perturbation would make 
such unstable orbit recede from r = rc 2 to r = ta, then reflect at r = ta, or move from 
r = tq 2 to r = re, then reflect at r = tb- The particle will move between r = ta and r = 
and will make a unstable choice on a movement direction due to the perturbation. Figure 6 
shows two cases numerically. 

II) When the energy of particle E equals the bottom value Eq\ of the effective potential, 
the particle can orbit on a stable circular orbit at r = r^. Or the particle orbits on a 
many- world bound orbit in the range ta < r < rp, where the ta and ro are the perihelion 
and aphelion distance, respectively. Figure 7 shows two cases numerically. 

C. Time-like escape geodesies 

When the particle energy is above the critical value (i.e. the peak value Ec 2 of the 
effective potential), the particle can orbit on a two-world escape orbit with a curly structure 
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FIG. 6: Examples of the unstable time-like circle orbit in the Bardeen space-time with Eq = 0.92, 
g = 0.7m, L = 3.5m and m = 1. 



1.00 




FIG. 7: Examples of the stable time-like circle orbit in the Bardeen space-time with E 2 = 0.90, 
g = 0.7m, L = 3.5m and m = 1. 



and cross the two horizons which is showed in Fig. 8a. When the energy of particle is much 
higher than the critical value, the escape orbit straightly deflects without curls, which is 
showed in Fig.8b. This means that the test particle coming from infinite would be reflected 
at a value of r and would not be able to reach r = 0, due to the infinite potential barrier 
at r = 0. on the other words the particle approaches the black hole from an asymptotically 
flat region, crosses the horizons twice and moves away into another asymptotically region. 



V. NULL GEODESICS 

For the null geodesic h = 0, we get the corresponding effective potential from Eq.flT] 

2mr 2 
[r 2 + g 2 )i' r v 

The behavior of the effective potential depends on the parameters g, L and the corresponding 
orbit equation is 

o ~o / *2tY17* 2 . Tj 2 , . 

r 2 = E 2 -(l- )-. (20) {eq : nullmo 

[r z + g z ) 2 r 
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FIG. 8: Examples of the two-world escape orbit in the Bardeen space-time with E 2 = 1.2, g = 0.7m, 
L = 3.5m and m = 1 (top) and the time-like escape orbit in the Bardeen space-time with E 2 = 7, 
g = 0.7m, L = 3.5m and m = l(bottom). 




FIG. 9: The behavior of the effective potential of the null geodesies for fixed L = 3.5m, m = l(left) 
and for fixed g = 0.8m, m = l(right). 

By using Eq.( l20|) and making the change of variable u~ l = r, we obtain the orbit equation 
for massive particle 

,du.~ E 2 2 2mu 
d<f> L 2 (l + u 2 g 2 )2 

By differentiating the Eq. (|2ip . we have 

d 2 u 3mu 2 . 

1^ +U= 7 T^J- 22 {eq = null or 

d(p 2 (l + u 2 g 2 )2 

We must solve the geodesic equations (12Tj) and (1221) numerically to investigate the null 
geodesies structure and how the space-time parameters influence on the null geodesies struc- 
tures in the Bardeen spacetime. We continue to follow the similar process of the Section 
IV. 
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A. Null bound geodesies 



From the effective potential curve for photons in Fig. 10, we can see that there are two 
different types of orbit when the energy E belows the peak energy value Ec- When the 
initial position is between ta and r#, the particle will move on a many- world bound orbit 
with the range of radius from ta to r#. When the particle initial position is on the right 
hand side of the potential barrier, the particle approaches rjj from an asymptotically flat 
region, then will be reflected to move away into another asymptotically region. These two 
kinds of orbits corresponding to the energy level are plotted on the right side of Fig. 10, 
respectively. 




FIG. 10: Examples of of the many-world null bound and escape geodesies in Bardeen spacetime 
with E 2 = 0.4, g = 0.6m, L = 3.5m and m = 1. 

B. Null circle geodesies 

When the energy E = Ec, The photon can orbit on a unstable circular orbit at r = tq, 
and the photo on such orbit will more likely recede from rc to ta crossing the horizons and 
will be reflected at ta or escape to the infinity on the other side of the potential barrier due 
the initial conditions and outside perturbation. Examples of such two kinds of orbits are 
shown in Fig. 11. 

C. Null escape geodesies 

When the energy E > Ec, The photon will be on the three different kinds of the escape 
geodesies which are shown in Fig. 12. We can see that when the energy level becomes 
larger from 0.6 to 7, the corresponding orbit changes from the two- world escape orbit to the 
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FIG. 11: Examples of of two kinds of unstable null circle geodesies in Bardeen spacetime with 
E 2 = 0.52, g = 0.6m, L = 3.5m and m = 1. 

escape orbit without intersection point, this means the particle approaches the black hole 
from an asymptotically flat region, crosses the horizons and then moves away into another 
asymptotically region. 



FIG. 12: Examples of three kinds of null escape geodesies in Bardeen spacetime with g = 0.6m, 
L = 3.5m, m = 1 and two energy levels E 2 = 0.6 and E 2 = 0.8. 



VI. CONCLUSIONS 



By analyzing the effective potential of massless and massive test particles in the Bardeen 
space-time, which describes a regular space-time and also represents a singularity free black 
hole for g 2 < (16/27)m 2 , where the parameter g represents the magnetic charge of the nonlin- 
ear self-gravitating monopole, and numerically simulating all possible orbits corresponding 
to all kinds of energy levels, we have found that there exist two kinds of bound orbits, one 
is close to the center of the black hole and crosses the two horizons, the other is outside 
the exterior horizon. The interesting result is that the planetary orbital precession direction 
is opposite and heir precession velocities are different, the inner bound orbit shifts along 
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counter-clockwise with higher velocity while the exterior bound orbit shifts along clockwise 
with low velocity, as shown in Fig. 5. We have also found two kinds of circular orbits, the 
inside one which closes to the exterior horizon is unstable and the outside one is a stable 
circular orbits, and two kinds of escape orbits. For the photon particle, there only exist one 
many-world bound orbit which can cross the inner- and out-horizons, one unstable circular 
orbit and three kinds of escape orbits. 
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